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Abstract 

In this work linear stability analysis of a magnetized dusty plasma 
with an anisotropic dust component having transversal motions much 
stronger than motions parallel to the external magnetic field, and 
isotropic light plasma components is described. Such a situation pre- 
sumably establishes in a shock compressed space dusty plasma down- 
stream the shock front. Oblique low-frequency magneto-hydrodynamic 
waves (u) <^ ojcdi ^cd being the dust cyclotron frequency) are shown to 
be undergone to the mirror instability. Consequences for nonthermal 
dust destruction behind shock fronts in the interstellar medium are 
discussed. 

PACS numbers: 32.80.Lg, 52.20.Hv, 52.30.-q, 52.35Py 



I. INTRODUCTION 

As shown in [H [2l [3l HI [5] the impurity of charged dust particles changes 
dynamical behavior of a plasma in the low-frequency limit where heavy dust 
particles dominate dispersion properties. Formally, the presence of an ad- 
ditional component increases the order of the dispersion equation and thus 
brings a new wave mode. One of interesting phenomena connected with the 
presence of dust in plasma is that due to a high inertia of dust grains such 
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complex plasma when brought once into an anisotropic state, can then be 
kept anisotropic much longer than an ion-electron plasma, and therefore can 
develop specific anisotropic instabilities. Such conditions are easily estab- 
lished in a magnetized dusty plasma subjected to shock waves. Indeed, after 
crossing a shock front with parallel magnetic field a dust particle starts gyrat- 
ing with the velocity V(i± = 3vs/4 relative to the light components (electrons 
and ions) (see discussion e.g. in [HI E]), where is the shock velocity. The 
longitudinal velocity component of dust particles remains however equal to 
the velocity dispersion of dust before the shock Vd\\ = Vdo Vs- As a rule, 
in these conditions the dust gyration period 2ii/ucd is much shorter than the 
dust drag time Td- oJcd^d ^ 1- For instance, in the interstellar dusty plasma 
behind a supernova shock with = 100 km s~^ 

30\Z\ 

UJcdTd ~ , (1) 

n 

where the drag was assumed to be due to direct ion-dust collisions as the 
Coulomb cross-section at these conditions is much smaller than the kinetic 
one, n is the plasma (ion) density, Z is the grain charge: ^ 3akT/e'^ - 
it can be as high as 10^ behind a shock with Vg = 100 km s~^; here Ucd = 
\Z\eB /mdC, r^^ = mpUadV / nid, nid, nip are the dust grain and the proton 
mass. Unless specified, we assume through the paper the characteristic size of 
dust grains of a = O.lyum, though our conclusions are weakly sensitive to the 
actual value of a since we consider the frequency range lower than the dust 
cyclotron frequency in the whole range of grain sizes. In particular, as Ucd'Td oc 
the condition UcdTd ^ 1 remains valid for the whole dust size spectrum 
ranging from a = 10 A to a = l/im (for the interstellar dust properties see 
the review [7]). It is therefore clear that the instability connected with the 
anisotropy of the dust component can develop succesfully on times t ^ oj~l 
provided the conditions for the instability are fulfilled. In this paper we shall 
describe an anisotropic instability in a plasma with isotropic electrons and 
ions of Tg ~ Tj, and anisotropic dust component with the transversal orbiting 
velocities much higher than the longitudinal velocities Vdi_ ^ Vd\\- 

The paper is organized in the following fashion. In Sec. II we derive 
the dispersion equation for magneto-hydrodynamic (MHD) waves in a three- 
fiuid plasma; in Sec. Ill we analyse solutions of the disperion relation and 
show that low-frequency waves are unstable; in Sec. IV consequences for dust 
destruction behind shock waves in the interstellar medium are discussed. 
Finally we close our paper with a short summary in Sec. V. 



2 



II. DISPERSION EQUATION 



Let us consider MHD waves in a plasma with cold dust particles gyrating 
with equal transversal velocities in an external magnetic field Bq = (0, 0, Bq), 
in which case their distribution function is /(v) oc S{vz)S{v± — vq), where vo 
is the transversal velocity of dust particles relative to the plasma component: 
Vq Vq = 3ws/4 for a plasma behind the shock wave moving with the velocity 
Vg] the wavevector is assumed to have the components k = (k_L,0, k||). Let 
us assume the phase velocity to fulfil the condition 

fT* < < VTe, (2) 

where vts and vti are thermal velocities of the electrons and ions respectively. 
In the low- frequency limit u <^ u^, where Ud is the ion cyclotron frequency, 
the equation of MHD waves is [8] 

here £,„,, e,,^ the corresponding components of the permeability. Un- 

der the validity of Eq. ([2]) and the condition that the ion gyration radius is 
smaller than the wavelength k±VTi <^ ood the components have the form 

^yy ^yy ^ ^yyi ^yy , ,2 ' 



u;2 ^2^2' 

where ojpe and Upi are the plasma frequencies of the electrons and ions, Uc 



is the cyclotron frequency of the electrons, Cs = \jTe/mi is the ion sound 
speed; the contribution of dust into the permeability components Syz and 
can be shown to be negligible. 

Let us assume the fiux of cold dust particles behind the shock front to be 
described by the distribution function 

/(v) = ^5(tg5K-^o), (5) 
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where n^o is the dust component density, Vq = 3fs/4 as defined above; the lon- 
gitudinal motions are thus explicitly considered to be negligible as compared 
to the shock propagation. This distribution function corresponds therefore to 
the extreme anisotropy, like if the temperatures ratio was Td±/Td\\ = oo for 
maxwellian dust distributions, and from general point of view in these con- 
ditions anisotropic instabilities are always expected to grow. The component 
Eyy has the form 

, .2 n=oo °? J „ 

4 = — E / ^^S— / dv^-'j'ni^Mv), (6) 



where 



F{v) 



VoVi_ 



\ k\\vA d6{v± - Vq) k\\v±d6{v^) 
1 ^ 6{v^) + 6{v± - fo) 



u J dv± u dvz 

^ (7) 

here x = k^v±/ucd, <^pd is the dust plasma frequency, Ucd, the dust cyclotron 
frequency. 

After simple algebra in Eq. (E]) one arrives at 



'yy 



E 



X 



Jn{x)J^{x) 



Let us consider low-frequency electromagnetic waves with u <^ Ucd, and the 
transversal wavelength significantly longer than the gyration radius of dust 
particles: k±v±/u}cd ^ 1- Then Eq. (|H]) can be rewritten in a simple form 



'yy 



^pd 

^cd 



k^^v? 



1 - 



2^ 



(9) 



The dispersion equation ([3]), with accounting (jlj) and ([9]) can be written 



as 



2^2 



fc^Q — k^c 



a;2 - klcl 



^2 7,2 



kj^ 



2^2 



where the notations c\ = BQ/Anpi, S = Pi/pd are introduced, pi, 
mass densities of the ions and dust. 



(10) 
are the 
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Equation fllOp is bicubic with respect to and can be rewritten as 



2{1 + S)n^ 



2 + (3 + 25/sin^ 6 + S{(3 + f) cos^ 6 



/3 + -W(l + sin2| 



2 



COS^ en'^ + -5/3/2 gij^2 Q ^Qg4 ^ _ ^^^^ 



Here 9 is the angle between the wavevector k and the external magnetic field 
Bo, n = uj/kcA, /3 = 2c2/4, / = v^/c\. 



III. RESULTS 



In the lack of dust Eq. ffTTj) describes two low-frequency electromagnetic 
modes present in a warm magnetized plasma: the fast and the slow magne- 
tosound waves. Indeed, when S = 0, Eq. f lTTj) degenerates and converges to 
a biquadratic one with the two solutions corresponding to the fast and slow 
magnetosound waves. When 6 and / = 0, Eq. ( ITT]) accounts the pres- 
ence of dust in the fast and slow waves in the limit u <^ Ucd'- the solutions 
are the fast and the slow waves loaded by the dust with the factor (1 + S). 
For gyrating dust particles (/ ^ 0) the order of the dispersion equation in- 
creases, and a third low-frequency mode emerges caused by the interaction 
of gyrating dust particles with the magnetic field of the wave. 

It is readily seen that this low-frequency mode is always aperiodically 
unstable, and corresponds to the mirror instability described in [9l [TOl [TT] . 
A typical value of the dust-to-gas mass ratio 6 in the interstellar medium 
is much lower than one (e.g. [7]), and as a rule f3 behind the shock front 
varies around one. At the same time / is of the order of the ratio of the 
shock velocity to the Alfven speed and can significantly exceed one. In the 
limit / — )■ oo the growth rate of the unstable low-frequency mode can be 
readily estimated asymptotically as F ~ f^^'^, while the frequences of the 
wave modes are 



4(1 + 5) 



sin^ 9 
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sm 



1 + S 



sin^ 29 



(3 



ni='- cos' 9. 



(12) 
(13) 



In the opposite limit / — > the solutions are 
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r ^ 



fV6 



sin 2^, 



(14) 
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for the growing mode, and 



1 



2 + /3 ± .\/(2 + /3)2 - 8/3(1 + 6) cos^ 6 +0{f), (15) 



4(1 + (5) 



for the fast and slow modes correspondingly. 

The effect of dust gyration on magnetosonic waves is illustrated on Fig- 
ures 1: Fig. la and lb show the frequency of the fast wave (a), and the 
growth rate of the third (low-frequency) mode (b) versus /; the slow mag- 
netosonic mode is rather unsensitive to /. Slowly gyrating and nongyrating 
dust particles (/ < 1) lead to a decrease of the phase velocity, which can 
be readely understood as due to a higher density of a dusty plasma than a 
dust-free plasma. 

As / grows the phase velocity increases, which is due to an enhance- 
ment of the transversal elasticity by the gyrating dust particles. Indeed, the 
transversal dust pressure can be defined as 



For strong shock waves, which are harmful for dust particles, Vq ^ ca, and 
the transversal pressure of dust component pdVQ becomes comparable (or 
exceeds) to the magnetic pressure Pic\/2 even though 6 = pd/pi <^ 1- 

Fig. 2 shows the dependence of the frequency Q and the growth rate F 
on plasma-/?. The phase velocity of the fast wave increases with /3 due to a 
joint action of the magnetic and thermal pressure. At small (3 the slow wave 
does not reveal significant deviations from the slow wave in an electron-ion 
plasma; deviations becomes considerable at sufficiently large /3. The growth 
rate of the unstable mode is rather unsensitive to the effects of the plasma 
thermal pressure. 

The mechanisms of the instability is similar to the mechanism of the 
mirror instability of the slow MHD wave in a plasma with anisotropic ion 
temperature Ti± > Ti\\. In our case, however, the isotropy of the electron 
and ion temperatures stabilize the slow and the fast waves, while on the 
other side the presence of a third cold gyrating dust component with an 
anisoptopic velocity distribution results in a new unstable mode. As expected 




(16) 
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the maximum of the growth rate T is reached when the angle between the 
wave vector and magnetic field is ir/A. 

IV. DISCUSSION 

We have shown in this paper that the presence of an anisotropic dust in a 
plasma results in development of the mirror instability in the low-frequency 
limit u ^ Ucd- The instability can be important for the fiows behind mag- 
netized shock waves in the interstellar medium: here dust particles keep 
significant transversal component during a long time Ucdt 3> 1 after crossing 
the shock front, and thus support development of the instability. The time 
of isotropization of dust grain velocities is of the order of Tj ~ ma/nmpadVo, 
where rup is the proton mass, ~ lO^^rrip is the grain mass, cr^ ~ 10~^° 
cm^, the grain cross-section, n is plasma density; this gives ~ s 
for typical value of the shock velocity Vq ~ 100 km s~^, thus resulting in 
^cdTi ~ 10 ''/ri ^1. It is worth noting that when direct collisions between 
protons and dust grains dominate, UcdTi increases when the grain radius a 
decreases, such that the dust component remains anisotropic on the growth 
time of the instability for the whole range of grain sizes. Grain-grain col- 
lisions remains unimportant in isotropization of the transversal dust veloc- 
ities: indeed, the characteristic time of grain-grain collisions is of the order 
Tgg ~ l/udadVod, such that Tgg ~ pTi/pd ~ IOOt/. Note that Coulomb colli- 
sions between charged dust particles are inefficient because the cross-section 
(jQ ~ Z^e'^ In A/ {nidv'^Y ~ 10~^^lnA cm^. In principle, fiuctuations of dust 
charge can enhance cxc (see, e.g. discussion in Ref. [HI [17]), however, they 
can be important only for dust grains with the charge Z ~ 1, while in the 
conditions with the charge Z ^ 100 what we are actually interested in, its 
fiuctuations cannot significantly enhance ac- 

As the normalized growth rate F depends only on 5, /3, / and the angle 6, 
the instability grows on a nearly Alfven crossing time ^ (0.3— 0.5) x A/c^ (see. 
Fig. lb and 2c). The upper limit of the wavelength for which isotropization 
can become important, and therefore the driving mechanism weakens and 
the growth rate drops, is determined from the condition F/cc^Tj < 1, and 
is of the order A ~ lO'^^CA/nvo cm. For the conditions typical behind the 
supernovae shocks in the interstellar medium ca ~ 30 km s~^, and the shock 
velocity Vq ~ 100 km s~^, this saturation wavelength is A ~ 3 x 10^^ cm. 
This value is approximately 3 times of the thicknes of the shell of a typical 
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supernova remnant when its expansion velocity is Vg = 100 km s~^. 

This means that the mirror instabihty induced by the anisotropy of dust 
velocities develops on all scales behind the shock front of a supernova rem- 
nant. As the growth rate decreases with the wavelength, the most fast grow- 
ing perturbations are expected to lie in the lower end of our approximation 
of the low-frequency limit, i.e. on several dust gyration radii. It is shown in 
Ref. [12] that on nonlinear stages the growing instability can form holes in 
the background magnetic field of several dust gyroradius in size, as observed 
in the experiment [T3] . 

From this point of view it can have great importance for dust destruc- 
tion behind strong shock waves in the interstellar medium. It is commonly 
known that strong {vg > 100 km s~^) radiative supernovae shocks efficiently 
destroy dust grains via sputtering (see e.g. p3| [T5]). The destruction ef- 
ficacy is mostly supported by betatron acceleration of the dust transversal 
motions / B{t) =const through the growing magnetic field in a radiatively 
cooling and compressing plasma: B cc p; p± is the transversal momentum of 
a grain. This means that the destruction mechanism works efficiently only 
behind the fronts with a planar magnetic field. In the fiows where the mir- 
ror instability operates and violates regular structure of magnetic field, the 
growing magnetic holes suppress acceleration of dust grains and therefore 
their destruction can be reduced. The exact value of the destruction rate in 
such conditions depends on many factors, such as the number of magnetic 
holes, their sizes and magnetic field strength, their lifetime, and can be only 
estimated in nonlinear numerical consideration. One should also mention 
that dust charge fiuctuations [161 [17] can result in acceleration/decceleration 
of dust particles, which in turn can affect the mirror instability. However, 
such effects apparently are important only when dust charge is small, close 
to Z ~ 1. In the conditions behind strong shock waves with highly charged 
dust grains \Z\ ~ 10^ — 10^ the increase of Eyy due to charge fluctuations, 
and therefore their influence on the mirror instability, seems to be small. 

IV. SUMMARY 

In this paper we have analyzed stability conditions of a dusty plasma with 
isotropic light components (ions and electrons) and anisotropic dust particles, 
with transversal velocities much higher than velocities along the external 
magnetic fleld. We have shown that 



8 



• oblique low-frequency {u ^ Ucd) magneto-hydrodynamic waves are 
unstable against the mirror instability; 

• this instability can be important in the interstellar dusty plasma behind 
strong shock waves from supernovae explosions, where the pressure of 
the gyrating dust particles can be comparable to the magnetic pressure 
behind the shock front. In these conditions the instability develops 
practically in the whole range of spatial scales, though the most fast 
growing wavelengths correspond to the lengths of several dust gyration 
radii; 

• the instability operates to destroy regular structure of the magnetic 
field behind the shock front, and can reduce the destruction efficiency 
of the interstellar dust by shock waves. 
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Figure 1: Q of the fast MHD wave (a) and F of the growing low- frequency 
mode (b) versus the gyration rate f'-S = 0.02, 0.05, 0.1 from bottom to top; 
13 = 0.01, cos^^ = 0.5. 
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(c) 



1 2 3 4 5 



Figure 2: The frequency Q of the fast wave (a) and the slow wave (b), and 
the growth rate F of the unstable mode (c) versus (3: 6 = 0, 0.02, 0.05, 0.1 
from bottom to top; k = 10, cos^6' = 0.5; as F = at 5 = curve 1 is not 
shown on pannel (c). 
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